function yp=cjxxfun(t,y,flag,q,m)
yp=[y(4);
y(5);
y(6);
q*(10*y(3)+1)*y(5)/m-q*((-5*sqrt(y (1).^2+y(2).^2))*cos(pi/4))*y(6)/m;
q*((-5*sqrt(y(1).^2+y(2).^2))*sin(pi/4))*y(6)/m;
-q*((-5*sqrt(y(1).^2+y(2).^2))*sin(pi/4))*y(5)/m];
end